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Introduction 

A correction to the previously published article “New Form of Kane’s Equations of Motion for 
Constrained Systems” 1 is presented. Misuse of the transformation matrix between time rates of 
change of the generalized coordinates and generalized speeds (sometimes called motion variables) 
resulted in a false conclusion concerning the symmetry of the generalized inertia matrix. The 
generalized inertia matrix (sometimes referred to as the mass matrix) is in fact symmetric and 
usually positive definite when one forms nonminimal Kane’s equations for holonomic or simple 
nonholonomic systems, systems subject to nonlinear nonholonomic constraints, and holonomic 
or simple nonholonomic systems subject to impulsive constraints according to Refs. 1, 2, and 
3, respectively. The mass matrix is of course symmetric when one forms minimal equations for 
holonomic or simple nonholonomic systems using Kane’s method as set forth in Ref. 4. 
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Symmetry of the Generalized Inertia Matrix 

Let S be a simple nonholonomic dynamical system, and let 7Z be an inertial frame of reference 
in which the configuration of S is described by a set of n generalized coordinates qi, . . . , q n . The 
velocity of a generic particle P of the system relative to 1Z can be written as 

n 

U V P = n ^r (?, t)u r + *vf (?, t) (1) 

r= 1 

where the generalized speeds ui,..., u n are scalar variables satisfying some nonholonomic con- 
straint relations linear in u\, . . . , u n , and 7 ' v(\ . . . ff v 7> are the corresponding holonomic partial 
velocities of the system. The generalized speeds also satisfy the kinematical differential equations 


Q = C{q,t)u + D(q,t) 


( 2 ) 


The system’s kinetic energy relative to 7 Z is given by 1 

\ 

K = - u T M (q, t)u + N(q, t)u + R(q, t) (3) 

where M G M nxn is symmetric, N e M lxn , and R e M. In most cases M is positive definite; 
it can be nonnegative definite (positive semidefinite) when, for example, one neglects a central 
principal moment of inertia of a rigid body belonging to S, such as a slender rod, or when one uses 
a particle with no mass to represent a point whose motion must be known. The rth nonholonomic 
generalized inertia force F* can be written in terms of the holonomic generalized inertia forces F* 
and in terms of K according to Eqs. (4.11.4) and (5.6.6) respectively of Ref. 4, 


n—p 


p* _ jp* 


E F U- A - 


5=1 

n r 


E 

5=1 


d / dK \ _ OK 
dt \ dq s ) dq s 


n—p 


W S r + E ^^SjP+k-^kr 


(■ r = l,...,p) 


(4) 


k = 1 


where the elements of W have the same meanings as in Eqs. (2.14.5) of Ref. 4, 


q = Wu + X (5) 

The matrix W is thus in fact identical to C(q, t); it is wrongly taken to be C~ 1 (q, t ) following Eq. 
(27) in Ref. 1, leading to errors in Eqs. (31), (32), (35), and (36), and to the incorrect conclusion 
that the generalized inertia matrix Q can be asymmetric. A rectification of the errors is provided 
in what follows. 
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The observation that generalized inertia forces are linear in u is made in Ref. 1 and will be 
revisited shortly. Consequently F*, a column matrix whose elements are F*, can be written in the 
form 

F* — —Q(q, t)u — L(q, u, t) ( 6 ) 

In Ref. 1 the nature of matrix Q is considered by using the kinetic energy of the system. Upon 
appealing to Eqs. (3) and (5) it can be seen that 


d K d K du r 7-1 

-7T7- = T7T = [u t M + N]W 
oq au oq L J 


Therefore, defining the matrix A 2 as in Eq. (16) of Ref. 1 


(7) 


A 2 


I 


( 8 ) 


yields the following matrix representation of Eqs. (4) 

A 2 F * = -A 2 W t [W~ t Mu + W~ t N t ] - K q T ^j 

= -A 2 W t (w~ t Mu + j f [W~ t M] a - ^ [W~ t N t ] - K^j (9) 

where 

K,= [K„...K,J = ^ (10) 

If Eq. ( 6 ) is multiplied by A 2 and compared with Eq. (9), we obtain 

Q = W t W~ t M = M (11) 

L = W T ^ [W~ t M] u+j t [W~ t N t ] - K q T ^j (12) 

Equations (11) and (12) correct Eqs. (35) and (36) of Ref. 1. The foregoing analysis shows that 
the generalized inertia matrix Q is symmetric and, as noted previously, usually positive definite. 
In addition to its role in Ref. 1, Q plays a central part in forming nonminimal Kane’s equations for 
systems subject to the constraints treated in Refs. 2 and 3. 


Alternative Demonstration of Symmetry 

It is noted in Ref. 4 that frequently it is inefficient to use K to construct F* or F*; likewise, 
readers of Ref. 1 are made aware that the use of K to form 0 or L is undesirable. The definition 
of F* given in Eqs. (4.11.2) and (4.11.3) of Ref. 4 provides alternative means for showing that 0 
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is identical to M, and for constructing M and L. When S is made up of particles Pi, , P v , 


F* = - ^2 m i n v Pi • n a Pi (r = 1, . . . , n) (13) 

i — 1 


where m, is the mass of Pi, and where the acceleration of P, in 1Z, denoted by R a Pi , can be written 
in terms of partial velocities by differentiating Eq. (1) with respect to t in 1Z as shown in Refs. 5 
and 2 . 

n 

r = 1 
n 

= ^2 1Z vf i u r + n = (14) 

r= 1 


('ll 
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where the remainder term of the acceleration is defined as 


n h l , 

Tl n P A ^ ' d 


'll 


a, = 


E 

r — 1 
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— V 

(It 


n^Pi 


u ' + it v * 


(i = !,...,//) 


(15) 


Substitution from Eqs. (14) into (13) yields 


p = -E 




2=1 

n / ^ 


E 

. S=1 


n v s Pi ft s + p a Pi 


y y, rrq v r R ‘ • | m. 


— ^ ///, /x v/’ • ^a t Pi (r = 1, . . . , n) (16) 


S=1 \ 2=1 


2=1 


It is immediately clear that F* is linear in p , .... u n ; each time derivative of a generalized 
speed is multiplied by an inertia coefficient defined as 


V 


A V P V P / \ 

m rs = 7 v m ?: v r 1 • v s 1 = m sr (r, s = 1, . . . , n) 
2=1 


(17) 


A matrix M whose elements are m rs is referred to as a generalized inertia matrix or mass matrix. 
Taken together, Eqs. (16) and (17) show conclusively that M is symmetric by definition when 
obtained with Kane’s method as set forth in Ref. 4. Although the present demonstration involves a 
holonomic system, the result is readily shown to apply also to a simple nonholonomic system; the 
interested reader is referred to Eqs. (5.5.3) and (5.5.4), and Problem 10.11 in Ref. 4. The form of 
m rs is dependent on one’s choice of generalized speeds because the partial velocities are dependent 
on this choice; however, M is symmetric for any valid choice whatsoever. The matrix M is in fact 
identical to Q that appears in Eq. ( 6 ) and thus a symmetric mass matrix appears in the nonminimal 
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Kane’s equations presented in Refs. 1, 2, and 3. The property of symmetry does not necessarily 
extend to a mass matrix obtained with other variants of Kane’s method, such as the lower triangular 
matrix obtained with the method presented in Ref. 6. 

It is important to note that not all methods of deriving dynamical equations guarantee a sym- 
metric mass matrix. For example, it is well known that Newton-Euler methods do not guarantee 
symmetry. Although Lagrange’s equations always produce a symmetric mass matrix, the same 
cannot be said in general for their cousins, the Boltzmann-Hamel equations. 7 

As indicated on p. 150 of Ref. 4, M is the matrix in Eq. (3) associated with the portion K 2 of 
kinetic energy that is quadratic in the generalized speeds 


K 2 


1 T 

-u T Mu 


1 

2 


n n 

EE TYl rs U r U s 

r = 1 s = 1 


(18) 


A comparison of Eqs. (16) and (6) reveals that the elements of matrix L are given by L r = 
Y^i = i m i U v r Pi ' Rsi t' ( r — 1 ,n). As pointed out in Ref. 1, it is advisable to construct F* 
according to the procedure set forth in Sec. 4.1 1 of Ref. 4 rather than by forming the matrices M 
and L. 
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